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ABSTRACT 



[ Several short-period Jupiter-mass planets have been discovered around nearby 

QQ I solar-type stars. During the circularization of their orbits, the dissipation of tidal 

' disturbance by their host stars heats the interior and inflates the sizes of these planets. 

■ Based on a series of internal structure calculations for giant planets, we examine the 

■ physical processes which determine their luminosity-radius relation. In the gaseous 
envelope of these planets, efficient convection enforces a nearly adiabatic stratification. 
During their gravitational contraction, the planets' radii are determined, through the 
condition of a quasi-hydrostatic equilibrium, by their central pressure. In interiors 



o 

^ [ of mature, compact, distant planets, such as Jupiter, degeneracy pressure and the 

^ ' non-ideal equation of state determine their structure. But, in order for young or 



intensely heated gas giant planets to attain quasi-hydrostatic equilibria, with sizes 
comparable to or larger than two Jupiter radii, their interiors must have sufficiently 
high temperature and low density such that degeneracy effects are relatively weak. 
Consequently, the effective polytropic index monotonically increases whereas the 
central temperature increases and then decreases with the planets' size. These effects, 
along with a temperature-sensitive opacity for the radiative surface layers of giant 
planets, cause the power index of the luminosity's dependence on radius to decrease 
with increasing radius. For planets larger than twice Jupiter's radius, this index is 
sufficiently small that they become unstable to tidal inflation. We make comparisons 
between cases of uniform heating and cases in which the heating is concentrated in 
various locations within the giant planet. Based on these results we suggest that 
accurate measurement of the sizes of close-in young Jupiters can be used to probe their 
internal structure under the influence of tidal heating. 
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1. Introduction 

One of the surprising findings in the search for planetary systems around other stars is the 
discovery of extrasolar planets with periods down to 3 days (Mayor & Queloz 1995). Nearly all 
planets with period less than 7 days have nearly circular orbits. In contrast, known extrasolar 
planets with periods longer than 2-3 weeks, have nearly a uniform eccentricity distribution. The 
shortest-period planets and their host stars induce tidal perturbations on each other. When these 
disturbances are dissipated, angular momentum is exchanged between the planets and their host 
stars, leading toward both a spin synchronization and orbital circularization (Rasio et al. 1996). 

Bodenheimer, Lin, &: Mardling 2001 (hereafter Paper I), considered the effect of tidal 
dissipation of energy during the synchronization of these planets' spin and the circularization of 
their orbits. In that analysis, they compute a series of numerical models for the interior structures 
of weakly eccentric Jovian planets at constant orbital distances under the influence of interior tidal 
heating and stellar irradiation. In these previous calculations, the interior heating rate per unit 
mass was imposed to be constant in time and uniformly distributed within the planet. Under these 
assumptions, they showed that Jovian planets can be inflated to equilibrium sizes considerably 
larger than those deduced for gravitationally contracting and externally heated planets. For the 
transiting planet around HD 209458, they suggested that provided its dimensionless dissipation 
Q- value is comparable to that inferred for Jupiter (Yoder & Peale 1981), a small eccentricity 
(e ~ 0.03) would provide adequate tidal heating to inflate it to its observed size (Brown et al. 
2001). Since the orbital circulation time scale is expected to be shorter than the life span of the 
planet, they also suggested that this eccentricity may be excited by another planet with a longer 
orbital period. The prediction of a small eccentricity and the existence of another planet are 
consistent with existing data (Bodenheimer, Laughlin, Sz Lin 2003). 

There are at least two other scenarios for the unexpectedly large size of HD 209458b. Heating 
by stellar irradiation reduces the temperature gradient and the radiative flux in the outer layers of 
short-period planets. This process could significantly slow down the Kelvin-Helmholtz contraction 
of the planet and explain the large size (Burrows et al. 2000). However, even though the stellar 
flux onto the planet's surface is 5 orders of magnitude larger than that released by the gravitational 
contraction and cooling of its envelope, this heating effect alone increases the radius of the planet 
by about 10%, not by 40% as observed (Guillot & Showman 2002). 

An alternative source is the kinetic heating induced by the dissipation of the gas flow in 
the atmosphere which occurs because of the pressure gradient between the day and night sides 
(Guillot Sz Showman 2002). In order to account for the observed size of the planet, conversion of 
only 1% of the incident radiative flux may be needed, provided that the dissipation of induced 
kinetic energy into heat occurs at sufficiently deep layers (tens to 100 bars). Showman Sz Guillot 
2002 suggest that the Coriolis force associated with a synchronously spinning planet may induce 
the circulation to penetrate that far into the planet's interior, and that dissipation could occur 
through, for example, Kelvin-Helmholtz instability. A follow-up analysis suggests that this effect 
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may be limited (Burkert et al. 2003; Jones & Lin 2003). 

In order to distinguish between these three scenarios, the effect of tidal heating for planets 
with modest to large eccentricities was considered. In a follow-up paper (Gu et al. 2003, hereafter 
Paper II), we showed that the size (Rp) of compact Jupiter-mass planets slowly increases with 
the tidal dissipation rate. For computational simplicity, wc adopted a conventional equilibrium 
tidal model which describes the planets' continuous structural adjustment in order for them to 
maintain a state of quasi-hydrostatic equilibrium in the varying gravitational potential of their 
orbital companion. In this prescription, a phase lag into the response is introduced to represent 
the lag being proportional to the tidal forcing frequency and attributable to the viscosity of the 
body. The phase lag gives rise to a net tidal torque and dissipation of energy and the efficiency 
of the tidal dissipation can be parametrized, whatever its origin, by a specific dissipation function 
or Q-value (quality factor) (Goldreich &; Soter 1966). External perturbation can also induce 
dynamical tidal responses through the excitation of g modes (loannou Sz. Lindzen 1993a,b) or 
inertial waves (Ogilvie & Lin 2004) which can be damped by viscous dissipation in the interior 
(Goldreich & Nicholson 1977) or radiative and nonlinear dissipation in the atmosphere of the 
planet (Lubow et al. 1997). 

In both the prescription for equilibrium and dynamical responses, the tidal dissipation rate is 
a rapidly increasing function of the planets' radius. But, their surface luminosity increases even 
faster with Rp such that, planets with relatively small eccentricities and modest to long periods 
attain a state of thermal equilibrium in which the radiative loss on their surface is balanced by the 
tidal dissipation in their interior. For planets with short periods and modest to large eccentricities, 
the rate of interior heating is sufficiently large that their Rp may inflate to more than two Jupiter 
radii. In this limit, the surface luminosity of the planet becomes a less sensitive function of its 
Rp and the eccentricity damping rate is smaller than the expansion rate of the planet so that 
the increases in their surface cooling rate cannot keep pace with the enhanced dissipation rate 
due to their inflated sizes. These planets are expected to undergo runaway inflation and mass 
loss. We suggested that the absence of ultra-short-period Jupiter-mass planets with P < 3 days, 
which corresponds to an orbital semi-major axis a of 0.04 AU, may be due to mass loss through 
Roche-lobe overflow resulting from such a tidal inflation instability. Other scenarios have been 
proposed to explain the lack of Jupiter-mass planets with P < 3 days: truncation of inner part 
of a disk (Lin et al. 1996; Kuchner & Lecar 2002), orbital migration due to the spin-orbit tidal 
interaction between the close-in planet and the parent star (Rasio et al. 1996; Witte &: Savonije 
2002; Paetzold & Rauer 2002; Jiang et al. 2003), and Roche-overflowing planets with the help of 
disk-planet interaction but without tidal inflation (Trilling et al. 1998). 

In this contribution, we continue our investigation on the internal structure of tidally 
heated short-period planets. The main issues to be examined here are: 1) how does tidal energy 
dissipation actually lead to the expansion of the envelope? 2) how does the internal structure of 
the planet depend on the distribution of their internal tidal dissipation rate? and 3) what are the 
important physical effects which determine the tidal inflation stability of the planets? These issues 
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are important in determining the mass-radius relation of short-period planets which is directly 
observable. 

Structural adjustments may also modify the efficiency of dissipation and the planets' Q-value 
for both equilibrium and dynamical tides. In the extended convcctivc envelopes of gaseous giant 
planets and low-mass stars, turbulence can lead to dissipation of the motion that results from 
the continual adjustment of the equilibrium tide. However, the turbulent viscosity estimated 
from the mixing-length theory ought to be reduced by a frequency-dependent factor owing to the 
fact that the convective turnover time scale is usually much longer than the period of the tidal 
forcing. Based on the present-day structure of Jupiter, Goldreich &; Nicholson (1977) estimated 
Q 5 X 10-^^. However, within the intensely heated (by tidal dissipation) short-period extra solar 
planets, convection is expected to be more rigorous with higher frequencies whereas their tidal 
forcing frequencies are smaller than that of Jupiter. The Q-value for the equilibrium tide within 
extra solar planets is likely to much smaller than that within Jupiter. 

For the dynamical response of short-period extra solar planets planet, the forcing frequencies 
are typically comparable to their spin frequencies but are small compared to their dynamical 
frequencies. Convective regions of the planets support inertial waves, which possess a dense or 
continuous frequency spectrum in the absence of viscosity, while any radiative regions support 
generalized Hough waves (loannou &; Lindzen 1993a, b). Inertial waves provide a natural avenue 
for efficient tidal dissipation in most cases of interest. The resulting value of Q depends, in 
principle, in a highly erratic way on the forcing frequency. Since the planets' spin frequency 
adjusts with their Rp, which is a time-delayed function of the tidal dissipation rate within them, 
the efficiency of tidal dissipation may fluctuate while the overall evolution is determined by a 
frequency-average Q-value (Ogilvie &; Lin 2004). 

In §2, we briefly recapitulate the basic equations which determine the quasi-static evolution 
of the planets' structure. We show the simulation results for inflated giant planets in the case 
of constant internal heating per unit mass in §3 and analyze the results, in the Appendix, in 
terms of polytropic models which allow us to conclude that the onset of the tidal runaway 
inflation instability is regulated by a transition in the equation of state for the interior gas from a 
partially degenerate/non- ideal state toward a more ideal-gas state. We examine the dependence of 
planetary adjustment on different locations of the energy dissipation in §4. Finally, we summarize 
the results and discuss their observational implications in §5. 



2. The Planetary Structure Equations and Numerical Methods 

The internal structure of an inflated Jupiter is constructed in this paper with the same 
numerical scheme as that used in Paper I. The code employs the tabulated equation of state and 
adiabatic gradient Vad described by Saumon, Chabrier, &: Van Horn (1995) and the values of 
opacities are derived from those provided by Alexander &; Ferguson (1994). With this scheme. 
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we compute distributions of thermodynamical parameters as a function of radius for a spherically 

symmetric planet (i.e. 1-D calculation). The surface temperature of the planet is assumed to be 
maintained by the irradiation from a solar-type star located 0.04 AU away. The structure of the 
planet is assumed to be hydrostatic and not to be affected by rotation since the nearly synchronized 
short-period planets spin at least several times slower than Jupiter and their rotational energy is 
quite small compared to their gravitational energy (paper I) 

The models include a tidal heating function, whose physical basis is not well understood. 
Thus in §3 we consider simulations in which the planet has a constant internal heating which is 
uniform in mass. In §4 we consider a set of models in which the heating is non-uniform in mass. 
As discussed by Ogilvie & Lin (2004) , different mechanisms may operate in radiative or convective 
layers. One should, in principle allow for three avenues of tidal dissipation: viscous dissipation 
of the equilibrium tide, viscous dissipation of inertial waves, and emission of Hough waves in the 
radiative zone. Rotational effects can also affect the behavior of tidal dissipation. In view of the 
uncertainties, we do not model specific mechanisms but simply parameterize the heating rate. 

In our calculations, we assume that convection is so efficient that the temperature gradient 
behaves adiabatically within a convective zone (Guillot et al. 2003). In other words, we solve the 
following equations for the radius r, the density p, the pressure P, the temperature T, and the 
intrinsic luminosity L (e.g. see Kippenhahn h Weigert 1990) 

(1) 



dm A-rrr^p^ 

dP _ Gm 

dm Airr^ ' 

dT GmT 



(2) 

V, (3) 



dm Anr'^P 

3 kCP 

V = "^rad = rrz 77 if radiative; V = Vad if convective, (4) 

_ _ s_dp^ 

dm ^ dt p dt ^ 
where m is the Lagrangian mass coordinate, 5 (= —{din p/dlnT)p) is the coefficient of thermal 
expansion at constant pressure, cp is the specific heat at constant pressure, e is the heating rate 
per unit mass, and V is defined as d\nT/dlnP. 

The above structure equations are solved simultaneously with an 1-D implicit Lagrangian 
scheme which uses the mass coordinate m as an independent variable. In the case where the 
heating rate is uniformly distributed with mass (i.e. e =constant), the results of our numerical 
calculations usually show that the planet interior is largely convective. Consequently, without 
the aid of the energy equation (5), equations (1) -(4) imply that the adiabatic assumption gives 
rise to a unique radial stratification for a given planet's size and mass, regardless of how strong 
the internal heating rate is or whether the planet is in thermal equilibrium. Equation (5) then 
indicates how fast the planet expands or contracts due to thermal imbalance. These adjustments 
proceed through a series of quasi-hydrostatic and quasi-thermal equilibria. 
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The initial condition is obtained from calculations for the formation phases of planets of the 
appropriate masses, as described in Bodenheimer et al. 2000. Just after accretion ends, those 
calculations show that the radius is about 2i2j, which is the initial condition used here. This 
value is somewhat uncertain, and it changes quite rapidly during the earliest part of the planet's 
cooling phase, but the exact value makes little difference here, because most of the models reach a 
thermal equilibrium which is practically independent of the initial state. 

For the inner boundary, we consider models with or without cores. For the models with 
cores, we assume that they have a constant density Pcore = 5 g cm~^. Their temperature is also 
assumed to be the same as that of the envelope immediately outside the core. The luminosity 
generated in the cores is assumed to be negligible. In most of our models, the central temperature 
is above ~ 3 x 10^ K, and the heavy elements in the cores are likely to be soluble (Paper I). 
Since we have already shown that core-less models lead to more inflated planetary structure, 
the planetary radii determined for the core models represent a lower limit. There are some 
uncertainties in the equation of state at the temperatures and pressures of the central regions, 
where the hydrogen-helium gas is partially degenerate and non-ideal. We show below that the 
degree of degeneracy and the non-ideal effects are important in determining the stability against 
the tidal runaway instability. 

Near the planet's surface, the gas becomes radiative when < Vad, where Vr is computed 
using equation (4) and Vad is obtained from a tabulated equation of state. The depth of the 
radiative zone is sensitive to the opacity. The temperature of the surface layer is sufficiently 
low for grains to condense. We adopt the standard opacity table (Alexander Sz Ferguson 1994) 
which is computed for the interstellar gas. However, in Jupiter's atmosphere the sedimentation of 
grains leads to a much reduced opacity which modifies the structure of the radiative zone. For 
short-period planets, however, the side facing the host star is heated to 1, 500 K which is near 
or above the sublimation temperature of most grain species. In addition, a large scale circulation 
flow may also modify the composition and the heat transport process near the surface. It remains 
to be determined whether grain sedimentation occurs (Burkert et al. 2003). In order to take into 
account this possibility, we consider some models under the assumption of extreme grain depletion, 
i.e. without any grain opacity, so that the effective opacity Kd is reduced from the normal opacity 
Ko (see Table 2). 

At the surface of the planet, we include the irradiation from the star in the boundary 
condition for the temperature Tgurf- We neglect the difference between the day and night sides 
of the planet. For most of the models, we set the irradiation temperature Tq = 1500 K, which 
corresponds to a planet with a 3-day period around a solar-type star. For two models, we also 
considered somewhat lower temperatures. At the photosphere, the normal boundary conditions for 
pressure P = 2GMp/3KciRp and for the intrinsic luminosity JC = AiraR^iTg^^^ — Tq ) are adopted. 
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3. Constant internal heating per unit mass 

In this section, we shall show the numerical results for an internal heating rate that is 
constant in time and is uniform in mass; that is, e in equation (5) takes a constant value. The 
calculations are carried out for three different masses: 0.63Mj, IMj, and 3Mj. The time span for 
each calculation covers a few 10^ years. At the end of each computation, the planets have attained 
a thermal equilibrium. In this state, the integrated energy term e in equation (5) balances the 
intrinsic radiated luminosity. Normally when a planet contracts and cools, the time-dependent 
terms in equation (5) dominate, and a thermal balance is never reached. The same is true even 
if the planet is strongly irradiated by the star. However the present results show that an internal 
tidal energy dissipation rate on the order of 10^^ to lO^'^ erg per second, depending on the mass, 
is sufficient to stop the contraction and cooling at some radius in the range of a few Rj. 

In Figure 1, we show that there exists a radius- luminosity relation for a Jovian planet 
with a given mass Mp. The planet's intrinsic luminosity (i.e. excluding the luminosity due to 

stellar irradiation) emitted from the photosphere is denoted by C. Discrete points shown in the 
Figure are our computational results, which are fitted by three quadratic curves associated with 
three different masses. For a 0.63Mj core-less planet, the solid curve in the left panel can be 
approximated by, 

log(/:/102^erg/s) = -11.7039(log(i?p/i?j))2 + 12.8994 log(i?p/i?j) - 2.12471. (6) 

For the same mass planet with a core (solid curve in the right panel) , we can approximate 

log(>C/102^erg/s) = -8.78501(log(i?p/i?j))2 + 10.004 log (i?p/i?j) - 1.05358. (7) 

Similarly, the model without a core (dashed curve in the left panel) and with a core (dashed curve 
in the right panel) for a IMj planet can be approximated by 

log(£/10^^erg/s) = -7.7909(log(i?p/i?j))2 + 10.3338 log (i?p/i?j) - 1.20809, (8) 

log(i:/102^erg/s) = -6.74289(log(i?p/i?j))2 + 8.92021 log{Rp/Rj) - 0.573055 (9) 

respectively. Finally, 3Mj models without a core (dotted curve in the left panel) and with a core 
(dotted curve in the right panel) can be approximated by 

log(/:/10^^erg/s) = -8.10772(log(i?p/i?j))^ + 10.6893 log(i?p/i?j) - 0.29812, (10) 

log(>C/10^^erg/s) = -10.525{log{Rp/ Rj)f + 12.3851 \og{Rp/Rj) - 0.508571 (11) 
respectively. 

For large values of Rp, these three curves all become flattened as the planet is inflated; i.e. 

as the planet expands, the slope (7 = d\og C / d\og Rp) decreases. In the core-less cases, 7 = 5 at 
Rp = 2.175Rj for Mp = 0.63Mj, Rp = 2.199Rj for Mp = IMj, and Rp = 2.243Rj for Mp = 3Mj. 
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In the core cases, 7 = 5 at i?p = 1.927Rj for Mp = 0.63Mj, Rp = 1.953Rj for Mp = IMj, and 
Rp = 2.243i?j for Mp = 3Mj. Also, the luminosity for a given Rp is an increasing function of Mp. 

The results in Figure 1 do not correspond to the exact thermal equilibrium solutions. 
Although the planets quickly establish a hydrostatic equilibrium, C does not necessarily equal the 
internal heating rate. In these calculations, C is the instantaneous intrinsic luminosity associated 
with the corresponding Rp. We use various initial conditions for Rp and e for our computation, 
but the jC vs Rp curves remain nearly identical to those shown in Figure 1. This invariance is due 
to the adiabatic structure as we have stated in the previous section. The number marked on each 
data point is the degree of electron degeneracy at the center of the planet in the core-less cases or 
at the surface of the core in the core cases. The degree of electron degeneracy is expressed by the 
ratio of the Fermi energy to the thermal energy 

D = Epermi/kT ^ p/{6 X lO'^T^/^). (12) 

As shown in the figure, the degeneracy is lifted as the planet expands. 

In Figure 2, we plot the temperature against the pressure inside a planet with Mp = IMj. 
Each curve corresponds to a different value of Rp. While the upper-right parts of the curves, 
where T and P are higher, represent the conditions for planet interiors, the lower-left portions 
of the curves show the T and P distribution near the planet's surface. The thick line across the 
upper parts of the curves marks the "Plasma Phase Transition (PPT)" line which segregates 
hydrogen molecules and metallic hydrogen atoms. This line is approximately drawn based upon 
an extrapolation of a hydrogen phase diagram (Saumon et al. 1995). As indicated by the slope of 
the curves, Vq^ has a value ~ 1/3 in the interiors where hydrogen atoms are in metallic form, and 
then drops down to around 0.2 as we follow the curves to the outer region of the planet where the 
pressure decreases below 10 bars. This transition occurs near the base of the radiative envelope. 
The flattening of V in the radiative envelope results from the external heating due to the stellar 
irradiation. The overall behavior of the curves for core and for core-less cases is very similar. In 
the core-less case, a decrease in the central pressure Pc with the planet's size can be fitted with a 
following power law 

Pc (X (13) 

The above relation is less steep than what is expected from a purely polytropic structure in 
which Pc oc In addition, the central temperature Tg increases with Rp when Rp ^ 2.3Rj, 

but decreases with Rp when Rp ^ 2.3Rj. This variation of temperature is nearly independent of 
the strength of internal heating rates; that is, the variation of temperature profiles is almost the 
same for a given planet's size no matter whether the planet is expanding, contracting or even in 
a thermal equilibrium. This tendency arises because, as stated in the previous section, adiabatic 
profiles are imposed in the convection zone which amounts to 95% of the total radius. The change 
of internal temperature in this manner is not driven by a thermal imbalance, but it is caused by 
an intrinsic change in the physical properties of the fluid in the convection zone as a result of the 
reduction of electron degeneracy (see Appendix A. 2). 
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Figure 3 displays the mass density profiles as a function of the temperature for various planet 
sizes. The theoretical Plasma Phase Transition line is also marked in the plot. Resembling the 
flattened P — T curves in Figure 2, the steepening of the T — p curves for p less than about 10~^ 
g/cm^, which roughly corresponds to the region of radiative envelope, is caused by the stellar 
irradiation. We estimate the degree of degeneracy at the center of the planet with no core, and 
at the surface of the core for the planet with a core. As the planet with no core (with a core) 
expands from 1.777Rj to 2.826i?j {l.GRj to 3.9Rj), the value D dechnes from 21 to 9 (from 21 
to 5.3) due primarily to the central density decrease. We fit the relation between pc and Rp for 
the planet with no core by using the following power law 

Pc oc i?;^^ (14) 

As is the case for the Rp — Pc relation, which is described by equation (13), the power 1.6 in 
equation (14) is smaller than what is expected from a purely polytropic structure in which the 
power is 3. 

The interior model of Jupiter is usually approximated by an n = 1 polytrope, which results 
from the combination of Coulomb effects and electron degeneracy in the metallic hydrogen plasma 
(Hubbard 1984; Burrows ct al. 2001). This simplification is confirmed by the simulation as shown 
in Figure 4 which plots P/p'^ (in c.g.s units) as a function of the co-moving mass coordinate m 
for various planet sizes. The hypothetical plasma phase transition line is denoted by PPT and 
is marked by a thick line. The almost constant value of P/ p^ throughout the region of metallic 
hydrogen indicates a n 1 polytropic structure. As we follow these curves to the outer regions 
of the planet where the molecular hydrogen dominates, the slopes of these curves increase from 
zero, meaning that the polytropic index is increased. The fractional region of the n = 1 polytropic 
structure shrinks as the planet's size becomes larger. This tendency suggests that the interior 
structure at various evolutionary stages of an inflated planet does not proceed in a self-similar 
manner, but behaves differentially as a consequence of the evolution of the equation of state. 
This motivates us to employ the simple polytrope approach to sketch the physical properties 
of an inflated planet associated with different values of the polytropic index n. By fitting the 
curves in the region of convection zones in Figure 4 for the core-less cases with a polytropic 
index, we show in Table 1 that the polytropic index increases with Rp. As demonstrated by the 
polytropic approach compared to the simulations in Appendix, the increase in n, and therefore 
the decrease in the adiabatic exponent F (see equation (A19) and the right panel in Figure 10 in 
Appendix), with Rp for the interior structure arises from a migration of an equation of state from 
a more degenerate and non-ideal phase to a state with less degenerate and more ideal properties. 
Consequently, for a given planet mass, the rate of Rp change in response to the change of total 
entropy becomes more sensitive as n increases (see equation(A21) for the planet compressibility 
Cm and the description in Appendix); namely, a larger planet is more elastic with respect to 
the change of entropy and pressure. This explains the simulation results that the increase in the 
intrinsic luminosity JC is reduced as the planet expands as shown in Figure 1 and equation (A17), 
making a larger planet vulnerable to the tidal inflation instability if the dimensionless parameter 
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for the tidal dissipation Q remains unchanged during the inflation. 

Although the above analysis is based on the simulation for uniform heating rate per unit 
mass, the concept of planet elasticity depending on the planet's size can be extended to the 
off-center heating cases. This will be discussed in the next section. 

4. Dependence on the location of the tidal dissipation 

4.1. Prescription for tidal dissipation rate 

In general, the distribution of tidal energy dissipation is a function of the radius, i.e. e = e(r). 
As we indicated above, the actual functional form e(r) depends not only the location of wave 
excitation and propagation, but also on the nature of dissipation. In stars and gaseous planets, 
the tidal perturbation of a companion not only distort their equilibrium shape but also excite 
gravity and inertial waves (Zahn 1977, 1989; loannou & Lindzcn 1993a, b, 1994, Ogilvie & Lin 
2004). The equilibrium adjustment and the dynamical waves which propagate within the planets 
are dissipated either through convective turbulence in its interior or radiative damping near the 
surface (Zahn 1977, 1989, Goldreich & Nicholson 1977, Papaloizou &; Savonije 1985, Xiong et al. 
1997, Lubow et al. 1997; Terquem et al. 1998; Goodman & Oh 1997). Tidal dissipation of the 
equilibrium tide deep inside a convective giant planet might be small because the convective flow is 
so adiabatic that the eddy turnover time of convection is much larger than the tidal forcing period 
which is about 3 days for the close-in planets we are investigating, leading to the outcome that the 
dissipation is restricted to an area far from the center of the planet interior. However, the tidal 
dissipation might take place at greater rates than the conventional calculations estimate through 
resonance locking between the perturbing tidal potential and oscillation (Witte &; Savonije 2002), 
which might happen to the tidal dissipation inside a giant planet. Most importantly in the case of 
a convective hot Jupiter, the dissipation rate associated with inertial waves in resonance with the 
harmonics of tidal forcing might not be severely reduced in the regime of low viscosity because 
the dissipation is independent of the viscosity (Ogilvie & Lin 2004). The dissipation of the waves, 
which occurs throughout the planets' convective envelope, also leads to the deposition of angular 
momentum which in general leads to differential rotation (Goldreich &; Nicholson 1977, 1989). 
It is not clear how the interior of the gaseous planets and stars may readjust after the angular 
momentum deposition may have introduced a differential rotation in their interior and how the 
energy stored in the shear may be dissipated (Korycansky et al. 1991). 

In this paper, we are focused on the response of the planets' envelope to the release of energy 
associated with the tidal dissipation. To determine the actual dissipation distribution is beyond 
the scope of the present work. In light of these uncertainties, we adopt three ad hoc prescriptions 
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for the positional dependence in e(r) such that 

= (^^1 1 I 

where Mp is the planets' total mass, /3 is a power index, mo and Am determine the Gaussian 
profile. For /? > 0, ei corresponds to a concentrated dissipation near the planet's surface such as 
that due to the atmospheric nonlinear dissipation or radiative damping, £2 corresponds to intense 
dissipation within the planets' envelope such as that expected from turbulent dissipation and 
nonlinear shock of resonant inertial waves, whereas 63 is designed to high light the dissipation 
near some transition zone such as the convective radiative transition front, or near the corotation 
radius. In the notion of equilibrium tides with constant lag angle, the total tidal energy dissipation 
rate within the planet in its rest frame (Eggleton et al. 1998, Mardling k, Lin 2002, papers I and 
II) is 

eo = ^tide = - f ^1 (^)'\nlh,{e)-2nn,h,{e)+n'h,{e)\, (18) 



p 



where and Mp are the mass of the stars and their planets, ji = M^Mp/{M^, + Mp) 
is the reduced mass, ^p and Rp are the planets' spin frequency and size, n, a, 
e are the mean motion, semi major axis, and eccentricity of the planets' orbit, 
/i3(e) = (1 + 3e2 + 3eV8)(l - e'^)'^/'^, /14(e) = (1 + 15eV2 + 456^/8 + 5eVl6)(l - e^)-^ 
and /15(e) = (1 + 31eV2 + 255e4/8 + 185e6/16 + 25e8/64)(l - e^)'"'^/'^. In the limit of low 
eccentricities, equation(18) can be approximated as follows: 



e'^GM^Mp 



where 



— (t) (^) (^) 

and Qp is the dimensionless parameter which quantifies tidal deformation and dissipation of a 
planet (paper II). 

An equilibrium tidal model with a constant lag angle for all components of the tide may 
not necessarily be the most reliable model. All the uncertainties associated with the physical 
processes are contained in the Qp values and we shall parameterize our results in terms of it. 
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For a fiducial value, we note that the Qp-value inferred for Jupiter from lo's orbital evolution is 

5 X 10'' < Qp < 2 X 10'^ (Yoder & Peale 1981). With this Qp-value, orbits of planets with Mp and 
Rp comparable to those of Jupiter, and with a period less than a week, are circularized within the 
main-sequence life span of solar- type stars as observed. 

For identical forcing to spin frequency ratios f^, the magnitude of eo for dynamical tides 
has the same power law dependence on Rp as that for equilibrium tide in eq (18). But, the 
dissipation rate and the Q-value vary sensitively with fn which is modulated by the changes 
in Rp (Ogilvie & Lin 2004). Numerical calculation and analytic approximation show that the 
relevant frequency-averaged Q-value is comparable to that inferred for Jupiter and it may be 
asymptotically independent of the viscosity in the limit of small viscosity or equivalently Ekman 
number. 

In the present context, we assume eg is constant in time in most of our models. But, for 
models 18-20, we consider the possibility that the damping time scale for the eccentricity (r^) is 
longer than the thermal expansion time scale (tr) for the planets to inflate. The expression for tr 
is given by the equation (paper II) 

4 . ( QpMpa qp Mp Mq2Rj ciV'^ 

r« = ^r., where = ^^^j ^ 0.18 (^^ — — ^^^J . (21) 

Hence in these models the heating rate eo increases as R^ according to equation (18) with a 
constant eccentricity. 



4.2. Numerical models 

In order to explore the dominant effect of tidal inflation, we consider several models (see 
Table 2) for an IMj Jupiter being inflated from the same initial size of 1.9 Rj. The parameters 
for these models are chosen to represent a wide range of possibilities. 

In the first and second series of models, mild (strong) heating is deposited in different locations 
in models 1-8 (9-12). With models 13-15 in the third scries, we consider the possibility of smaller 
opacities, due to grains' sedimentation, in the radiative envelope. With a lower Tg, we consider 
the possibility of less intense stellar irradiation in models 16 and 17 in the fourth series. Simple 
forms of time- varying tidal dissipation are also considered in the fifth series of models 18-20. We 
consider the damping timescale of eccentricity is longer than the expansion time scale, so the 
heating rate in models 18—20 takes the tidal form e oc i?^. We only deal with the core- less cases for 
all of these models shown in Table 2 except for model 19 in which the simulations for both core 
and core-less models are carried out. 
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4-2.1. Dependences on the heating location 

In this series of models, we consider the effect of non-uniform heating on the planets' 
radius-luminosity relation. The relation which stays approximately invariant for different locations 
of heat deposit is the radius-adiabat relation (Stahler 1988). In Figure 5, we plot the density 
profiles of an IMj inflated giant planet in a thermal equilibrium. Model 1 &; 6 are represented by 
the dotted line, because their interior structure is almost the same for these two models as the 
heating is concentrated in the inner regions of the planet. Models 2 and 5 are represented by solid 
and dashed lines respectively. 

These models indicate that as the location of maximal heating moves from the planet's center 
toward the photosphere {ttiq/Mj = 0.05 — > 0.7 0.9999), the final equilibrium size of the planet 
decreases. This tendency indicates that less PdV work on expansion is achieved when the heating 
shell is closer to the photosphere because the heating is also more efficiently lost near the planet's 
surface. Equivalently, less specific entropy content is retained by the convective region of the 
planet. Given a mass, a size, and the same heating rate, a planet with the heating shell closer to 
its photosphere has larger intrinsic luminosity £ as a result of the difficulty in transporting entropy 
into the planet's interior via both radiative and convective transport. Our numerical results 
indicate that the region beneath the main heating zone continues to adjust slowly without gaining 
any entropy. Even without a substantial increase in entropy, we show in the next paragraph 
for models 11 and 12, that the inner convective region can swell in order to adjust to a new 
hydrostatic equilibrium. This adjustment is due to a drop in the boundary pressure in the limit 
that the spacial extent of the heating shell spreads significantly. Although this result seems to be 
in qualitative agreement with the radius-adiabat relation, it is still unclear that the relation should 
precisely hold for a giant planet heated in different locations because different equations of state 
and therefore different PdV work are involved for a given entropy input to different locations. 

Figure 6 displays the evolution of expansion for model 11 (upper panels) and model 12 (low 
panels). The three solid curves from right to left in each of the temperature and density diagrams 
are the profiles corresponding to the phases marked by the data points 1, 6, 11 in the Rp vs t 
diagram for model 11, and the phases labeled by the data points 1, 4, 10 in the Rp vs t diagram for 
model 12. The vertical dashed lines in the temperature and density diagrams denote the location 
of the center of the Gaussian heating zones: m/Mp = 0.9 for model 11, and m/Mp = 0.7 for model 
12. In model 11, a heating front is generated at m/Mp = 0.9 as is illustrated in the temperature 
diagram (i.e. the upper-middle panel) and then propagates inward. The front corresponds to the 
region with the temperature inversion. It also acts as a rarefaction structure as is shown in the 
density diagram. For comparison, we also draw the temperature and density profiles for model 9 
(dot-dashed curves) and for model 10 (dotted curves). In model 12, a heating front is generated 
at m/Mp = 0.7. As it travels inward (see the lower-middle panel), the front decreases the density 
analogous to a rarefaction wave (see the lower-right panel). 

It is interesting to point out that except for the data points from 1 to 4 in model 12, all of the 
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other data points in the Rp versus t diagrams correspond to a state of quasi-thermal equilibrium; 
i.e. the intrinsic luminosity C is equal to the dissipation rate. However the planet is not in a global 
thermal equilibrium because of the presence of a temperature inversion in the region just below 
the shell. This heat front very gradually propagates inward, and is associated with a gradual 
overall expansion. After quasi-thermal equilibrium is reached, the rate of expansion is very slow 
compared to its rate before this time (compare the data before and after point 4 in the lower left 
panel of Figure 6) . Clearly, an analogous expansion phase does not occur in the uniformly heated 
case (model 9). 

The underlying physical cause for expansion of Rp without any gain of internal entropy (or 
very little gain of entropy as a result of the tail of a Gaussian profile) can be described as follows. 
Prior to reaching a quasi-thermal equilibrium, a giant planet experiences a fast inflation due to a 
high dissipation rate focused in a spherical envelope (such as the stage before the data point 4 for 
model 12 in Figure 6). In models 11 and 12, about 10% to 30% of planet gas by mass lies above 
the heating zone. This gas is inflated significantly due to a strong underlying shell-heating source 
as prescribed by the model. The huge inflation on the top of the heating shell reduces the pressure 
and therefore causes the whole region beneath the heating zone to expand and to adjust to a new 
hydrostatic equilibrium. The expansion induced by a shell-heating source is less efficient than that 
induced by a uniform heating throughout the envelope since more energy is radiated away for a 
given Rp, as we have already stated in the previous paragraph. The surface radiation becomes 
more intensified as the heating zone approaches to the photosphere. 

After the entire shell-heated planet has reached its quasi-thermal equilibrium, the region 
around the heating front is not yet in the thermal equilibrium locally. This arises because the 
region just below the heating front (i.e. at the bottom of the temperature inversion region) gains 
entropy via convection from the bottom and through radiation from its top. But the region just 
above the heating front (i.e. at the top of the temperature inversion region) loses entropy due to 
radiation from its bottom and convection from its top. On the average, the net increase in entropy 
vanishes. Thus, the rarefaction heating wave is driven by a local entropy transport which averages 
to a zero net flux. Since the heating front is a result of the temperature inversion which transports 
entropy inward by radiative diffusion and causes the entropy disturbance, the time scale for the 
heating front to cross the planet's interior is roughly equal to 




where the radiative diffusivity I'rad = 'iacT'^ /Snp^Cp, and Ar is the width of the heating front. The 
reference values for Ar, T, k, and p shown in the above estimate are taken from the simulation 
for model 12. The large value of traref action is consistent with our numerical results that the 
shell-heating planet expands very slowly after it has reached the quasi-thermal equilibrium (see 
the upper- left panel and the stage after the data point 4 in the lower-left panel in Figure 6). In 
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reality the slow expansion of the planet due to the inward propagation of the rarefaction wave is 
unimportant, because the damping time for tidal dissipation is much short than tj-are factions ^'^id 
also trarefaction is much longer than the lifetime of the star. 

The C — Rp relation for model 12 before the planet reaches a quasi-thermal equilibrium can 
be fitted with a quadratic approximation: 

log(>C/102^ erg/s) = -9.36888(log Rp/Rjf + 11.9403(log Rp/Rj) - 1.33544, (23) 

which indicates that 7 decreases with Rp and 7 = 5 when Rp = 2.346i?j. In contrast to the 
uniform heating cases, in which the planet's interior gains entropy and expands, the degree of 
degeneracy at the planet's center, in model 12, remains essentially constant [D 17) during the 
evolution (see Figure 6). This tendency arises largely because the planet interior expands without 
gaining much entropy. This result implies that the flattening of C with Rp in the shell-heating 
cases is not due to the reduction of electron degeneracy of the metallic interior, but is primarily 
caused by the gradual phase transformation from non-ideal properties to the ideal state within 
and above the shell-heating zone. This interpretation is well illustrated in Figure 7 in terms of J5 
(= 1/r, see equation (A19) and the description in Appendix) for the data point 1 (solid curve) 
and the data point 4 (dotted curve): while in the interior J 6 hardly changes with J 5 indeed 
increases with Rp around and above the shell-heating zone. 



4-2.2. Variations in other model parameters 

In addition to the cases in which most of the dissipation is concentrated in a narrow spherical 
shell, we also considered the cases with relatively flat heating profiles for 69. In Model 7 and 8, we 
choose the power index to be /3 = 2. We found that the interior structure is almost the same for 
Models 7, 8, and 1. This similarity arises because the heating per unit volume peg concentrated in 
the high-density inner region of the planet in all these models, including Model 7 with (3 = 2. For 
much larger values of /3, the internal structure of Model 7 and 8 should resemble that of Model 2 
and 6 respectively. 

Figure 8 depicts the density profiles of the planet in a thermal equilibrium with the regular 
and the reduced opacities for grains. The solid and dashed curves represent Models 1 and 4 
respectively. In both models, = no- The corresponding Models 13 and 14 for k,^ = 10~^fi:o 
are marked with solid and dashed curves respectively. This reduction of opacity by a factor of 
10^ is implemented in the grain region where T < 2100 K. This region roughly corresponds to 
the radiative envelope of the planet. Since the radiative envelope is a bottleneck for the outward 
energy transport, the reduced opacity in the radiative envelope allows photons to escape from the 
planet much more easily. Consequently, the planet attains a much smaller size when a thermal 
equilibrium is established. The dashed curves have slightly smaller Rp than the solid curves as 
a result of more radiative loss for the shell-heating cases as reasoned in the previous paragraph. 
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Since the location of the heating shell is close to the planet's center, the simulation for Model 15 
is almost the same as that for Model 13. 

We also computed the thermal equilibrium solution for models with different irradiation 
temperatures Tg. In Figure 9, we plot the temperature and opacity distributions for Models 1, 
16, and 17 with solid, dashed, and dotted curves respectively. Since the interior structures for 
these three cases are almost the same, we only show the temperature and opacity profiles in the 
radiative envelopes, which axe affected by the stellar irradiation. The planets with less irradiation 
(Models 16 & 17) have slightly larger equilibrium sizes than the one with more irradiation (Model 
1). A direct comparison between Models 16 and 17 indicates that higher surface temperatures 
result not only in a slight increase in the density scale height of the planet's atmosphere, but also 
an enhanced opacity. Both effects cause Model 16 to attain a larger Rp than Model 17. 

But Rp for Model 1 is smaller than that for Model 16 even though the photospheric 
temperature for the former is higher than that of the latter. This difference arises because the 
temperature of the envelope in Model 1 is sufficiently high for most of the silicate grains to 
sublimate. Consequently, opacity in the radiative envelope in Model 1 is well below that in Model 
16. The correlation between Rp and k is already established in the previous discussion on Models 
13-15. However, the difference in the size between Model 1 and Model 16 is only about 1%, in 
accord with the usual notion that the size of an optically-thick object should not be strongly 
altered by external irradiation. 



4-2.3. Self consistent tidal heating models and inflation instability 

We also consider, in Model 18, a self consistent calculation in which the energy dissipation is 
proportional to Rp in accordance with equation (18). In this model, we assume that the energy 
dissipation is uniformly distributed in mass. The dissipation is normalized to eo = 5 x 10~^Lq at 2 
Jupiter radii (this rate corresponds to a = IORq and e ^ 0.18 in equation (19). Note that ~ 
in this case according to equation(21)). The final equilibrium size Rp = 1.857 Rj is slightly smaller 
than that {Rp = 1.977 Rj) for Model 1. The planet would be slightly bigger {Rp = 1.983 Rj) than 
that for Model 1 if the normalization constant were increased by 5%. Since 7 < 5 as i?p ;> 2Rj, 
the planet is expected to be vulnerable to the tidal inflation instability once its size expands 
beyond 2Rj. Our numerical results show that a further increase in the normalization constant by 
another 5% (therefore e ~ 0.2 and tr ~ 0.9Td when Rp = 2Rj and a = IORq by equation(21)) 
causes Rp to increase to at least 4 Rj. 

In comparison to the uniform heating in Model 18, we consider in Models 19 and 20 the 
possibility of non-uniform heating with m = 0.9 and m = 0.9999, respectively. Similar to Model 

18, the total heating rate is proportional to Rp. In Model 19 with the heating rate normalized 
to 2.1 times larger than 5 x IO^^Lq at Rp = 2Rj (this corresponds to e ~ 0.26 at a = IORq 
according to equation (19)), the planet is thermally unstable and swells from 2Rj to 3Rj in 
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about 2 Myrs (shorter than the eccentricity damping time scale 2.7 Myrs according to 

equation (20) evaluated at Rp = 3Rj). When the heating rate is normalized to 4 times larger than 
5 X IO^^Lq at Rp = 2Rj (this corresponds to e ~ 0.07 at a = 0.03 AU according to equation(19)) 
in Model 19, the planet is thermally unstable and expands from 2Rj to 3Rj in about 0.27 Myrs 
(the eccentricity damping time scale 0.15 Myrs and 1.17 Myrs evaluated at Rp = 3Rj and 
Rp = 2Rj, respectively). 

When the normalized heating is set to be 10~^ Lq for a = 0.03 AU and Rp = 2Rj, equation 
(19) is no longer a fair approximation and equation (27) in paper II gives e 0.294 for this heating 
rate. In this case the planet with or without a core expands from 2 Rj to 3.6 Rj in less than 
40000 years which is shorter than the eccentricity damping time ~ 47800 years for Rp = 3.6 Rj 
according to equation (12) in paper II. Since the planet's size 3.6 Rj is actually beyond the Roche 
radius Rl ~ 3 Rj for a = 0.03 AU and e = 0.294, the inflated planet would overflow the inner 
Lagrangian point in this case. 

The expansion rate is drastically reduced for Model 20 in which the dissipation is largely 
deposited at m = 0.9999, roughly the location of the radiation-convection boundary. A one-Jupiter 
mass planet with the normalized heating 5 x W^Lq at Rp = 2Rj can only reach the final size 
of 1.9 Rj. a normalized heating rate which is twice as big as 5 x 10~^Lq results in a final 
size ~ 1.95 Rj. The planet with the normalized heating four times as large as 5 x lO^^L© in 
Model 20 expands very slowly to 2 Rj from 1.9 Rj over several tens of million years, which is 
comparable with the eccentricity damping time scale 20 Myrs at a = IORq and is much longer 
than the eccentricity damping time scale (2± 1.17 Myrs) at a = 0.03 AU. The planet in Model 20 
obviously requires a larger normalized heating rate than Model 19 to reach the critical size beyond 
which the planet is thermally unstable in response to the R^ heating rate. 

5. Summary and discussion 

In this paper, we continue our investigation on the adjustment of a planetary interior as 
a consequence of intense tidal heating. As a giant planet's interior is heated and inflated, we 
showed in Paper II that its interior remains mostly convective. Efficient energy transport leads to 
an adiabatic stratification. With a constant heating rate per unit mass, we deduced an unique 
luminosity-radius relation regardless of how intense the heating rate is. The planet's luminosity 
increases with its radius. But the growth rate of £ is a decreasing function of Rp. At the same 
time, the tidal dissipation heats the interior of the planet at a rate which increases rapidly with 
Rp. At around 2 Rj, C can no longer sustain sufficient growth to maintain a thermal equilibrium 
with the tidal dissipation rate. Thereafter, the planet's inflation become unstable and it overflows 
its Roche radius and become tidally disrupted. 

Here we show that the change of luminosity during the planet's expansion is directly linked to 
the evolution of its interior, in particular, the equation of state. We employ the polytrope approach 
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to investigate the interior structure of an inflated giant planet. According to the simulation, 
interior profiles deviate away from P oc as the planet expands. The central temperature Tc 
increases and then decreases with the size of the planet. Also P^. and Pc are less steep functions 
of Rp than the polytrope theory with a constant polytropic index n indicates. All of these effects 
suggest that the planet interior does not evolve in a self-similar manner, but n gets larger as Rp 
increases. In conjunction with the numerical results that the degeneracy D decreases, and that Tc 
rises and then drops during the course of inflation, the result of a positive value of dn/dRp can be 
interpreted as a manifestation of a reduction in degeneracy during the expansion. 

We reason that the coefficient of thermal expansion S increases in response to a decrease in 
degeneracy and non-ideal effects, leading to an increase in n through equation (A19). Consequently 
the planet compressibility at constant mass Cm increases with Rp (see equation A21). This 
pattern can be translated into the phenomenon of a decrease in luminosity growth during the 
inflation, as a consequence of an one-to-one relation between K and JC in the case of uniform 
heating in mass under the conditions of the polytropic interior and hydrostatic equilibrium. We 
also compare the results between a planet with a core and without a core. To be inflated to the 
same size, a planet with a core, therefore possessing a larger gravitational binding energy, needs a 
larger intrinsic luminosity C than a planet with no core and the same mass. We also show that 
the opacity in the radiative envelope has a drastic effect on the final equilibrium size of an infiated 
planet: the size would be much smaller if grains are depleted in the radiative envelope. 

We also consider the possibility of localized tidal dissipation. Such a process may occur in 
differentially rotating planets or near the interface between the convective and radiative zones 
where the wavelength associated with dynamical tidal response is comparable to the density scale 
height. Localized dissipation may also occur through the dissipation of resonant inertial waves or 
radiative damping in the atmosphere. In the strong shell-heating models, the one-to-one relation 
between K and C disappears because of the existence of a radiative region caused by temperature 
inversion beneath the shell-heating zone. The unheated planet's interior in such cases might still 
be infiated due to a significant expansion of the gas above the heating zone, although the overall 
expansion rate is less efficient than that in the uniform heating cases as a result of a greater 
amount of radiative loss from the planet's photosphere. Without gaining entropy, the expanding 
interior cannot lift its degeneracy and therefore cannot increase its elasticity. However, the gas 
above the shcU-hcating zone can lift its non-ideal properties and hence enhance its elasticity, 
leading to a decrease in F in that region and thereby diminishing the growth rate of C as the 
planet expands. 

Finally, we consider the self consistent response, taking into account the modification of 
heating rate due to a planet's expansion. In this paper, we adopt a constant-Q prescription for 
equilibrium tides in which the tidal dissipation rate is assumed to be proportional to Rp. The 
results for the uniform heating model suggest that a young gaseous planet of 1 Mj without a 
solid core can be thermally unstable and inflated from 2 Rj to a size beyond 4 Rj if e ~ 0.2 at 
a = 10 Rq. If the dissipation rate is proportional to Rp, and if most of the tidal perturbation 
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is deposited at m/Mp = 0.9, a core-less young planet of IMj would be thermally unstable and 
inflated from 2Rj to a size beyond 3Rj for e > 0.07 at a = 0.03 AU or e > 0.26 at a = lOi?0. 

With the same heating concentration m/Mp = 0.9 and i?p-dependence in the dissipation rate, 
a young planet with a core at a = 0.03 AU with an initial eccentricity e > 0.294 can be inflated 
from 2 Rj to a size beyond its Roche radius. Wc have assumed that the convectivc flow still 
behaves adiabatically even though the heating shell causes a narrow radiative zone. However, 
the condition away from adiabaticity implies that the internal heat is not transported away as 
efficiently as in the case of adiabatic convection, leading to a more severe reduction in non-ideal 
properties of the gas and therefore an even faster decrease of 7 as the planet expands. 

Note that the Eddington approximation for the surface boundary condition is used in these 
models rather than more detailed frequency-dependent model atmospheres. This approximation 
is not necessarily valid for the strongly irradiated atmospheres studied here (Guillot & Showman 
2002). However it is unlikely to make much difference for the main results discussed here, namely 
the behavior of the planet's radius as a function of tidal dissipation energy. It could, however, lead 
to errors in other kinds of predictions, such as the radius as a function of opacity. 

The equilibrium tidal dissipation formula is based on an ad hoc assumption of a constant 
lag angle. In reality, the dynamical tidal response of a planet through both gravity and inertial 
waves near the planet's surface and convective envelope may be much more intense, especially 
through global normal modes. Their dissipation may provide the dominant angular momentum 
transfer mechanism for the orbital evolution and heating sources for the internal structure of 
close-in extrasolar planets. In the limit of small viscosity, the intensity of tidal dissipation is 
highly frequency dependent (Ogilvie & Lin 2004). When the forcing and response frequencies are 
in resonance, the energy dissipation rate is intense whereas between the resonances it is negligible. 
As the planets undergo structure adjustments, their spin frequency. Brunt- Vaisala frequency 
distribution, the adiabatic index, and equation of state also evolve. Since all of these physical 
effects contribute to the planets' dynamical response to the tidal perturbation from their host 
stars, their response and resonant frequencies are continually modified. The results in this paper 
indicate that the structure of the planet adjusts on a radiation transfer time scale which generally 
differs from the time scale for a planet to evolve through the non resonant region. In addition, the 
tidal forcing frequency also changes as the planets evolve toward a state of synchronous spins and 
circular orbits. Therefore, it is more appropriate to consider a frequency averaged tidal dissipation 
rate. In the limit of small viscosity, the frequency averaged dissipation rate converges (Ogilvie &; 
Lin 2004) such that the equilibrium tidal dissipation formula may be a reasonable approximation. 
Nevertheless, we cannot yet rule out the possibility that some close-in planets may attain some 
non resonant configuration and stall their orbital evolution. Therefore, accurate measurement of 
the sizes of close-in young Jupiters via planet transit surveys can be used to constrain the theories 
of tidal dissipation and hence internal structure for these objects. 
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A. A Polytrope Model 

In order to identify the dominant physical effects which determine the internal structure 
of planets, we consider an inflated giant planet that consists of a polytropic interior and a thin 
envelope. With this model we can construct analytic solutions. For this analysis, we assume that 
the polytropic interior comprises almost all of the mass and radius, and that the thin envelope 
overlying the interior is radiative and is composed of a non-degenerate ideal gas. Therefore, 
the polytropic equation P = Kp^^^^'"-, together with the condition for hydrostatic equilibrium, 
specifies the value of K for a given planet's mass Mp and a given planet's size Rp (Cox & Giuli 
1968): 

K = A;„Mpi-V"i?-i+3/", (Al) 

where kn is a function of n. All of thermodynamic quantities must be continuous across the 
boundary between the convective interior and the radiative envelope, such that 

n = ^ (A2) 

n = i^pJ+^Z" (A3) 

polytrope rad ^q^^^q j,a > V > 

where the subscript h denotes the values evaluated at the boundary. The magnitude of V polytrope 
is determined by the equation of state of the gas in the interior such that it is a function of n 
rather than directly dependent on the magnitude of Rp and Mp. Note that in equation (A4) we 
have assumed that C does not vary greatly across the radiative envelope after it emerges from the 
polytropic interior. This assumption should be a reasonable approximation for the thin radiative 
envelope so long as there is no localized source of intense heating there. The above four equations 
give rise to the relation 



PbK-b Kb \TZ, 

where TZ is the gas constant and jj, is the mean molecular weight. From the opacity table, we find 
that the opacity at the boundary may be roughly approximated by a power law 

«6 oc Tfc", (A6) 



with a = 4.08. 
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A.l. Completely Degenerate Interior 

The equation of state for a completely degenerate gas in the non-relativistic regime is given 
by P oc (f'l^ , which corresponds to the n = 2>/2 polytrope with a constant K. Therefore, equation 
(Al) leads to the well-known mass-radius relation for low- mass white dwarfs: 

Rp oc M-V^ (A7) 

which does not vary with C This independence is equivalent to the expression 

d In £ ... 



A.2. Partially Degenerate Interior with n = 1 Polytrope 

We now consider a model in which the planet has a sufficiently large Rp such its interior is 
partially degenerate with an n = 1 equation of state {P cc p^). Partial degeneracy occurs when 
D 1. Therefore by setting D = 1 in equation (12), the partially degenerate interior can be 
described by 

p = cT^/^, (A9) 

where c is a constant. Under these conditions, V ad = 1/3, which is consistent with the numerical 
solution for the region dominated by the pressure- ionized hydrogen atoms as shown in Figure 2. 

At the boundary, equations (Al), (A2), (A3), and (A9) uniquely determine the value of Tj, for 
a given planet's size Rp without having to consider the photosphere: 

With equations (A5) and (A6), it follows that 

C oc MpRl^-^^c^ '^ (All) 

Figure 3 shows that the degree of degeneracy at the center of the planet decreases from 21 to 
9 as its radius expands from 1.777 Rj to 2.826 Rj. The value c decreases as the degree of partial 
degeneracy is reduced (c 6 x 10~^ g cm~^ K~^/^ when Fermi energy equals thermal energy). 
Therefore, in accordance with equation (A9), the temperature of the planet can increase with, 
even though p decreases with, Rp. ^ Temperature indeed increases until Rp is increased up to 
2.3i?j. The temperature then decreases with Rp after that, as seen in our numerical results. 



^Extrapolating from equation (A9) to the case of high degree of partial degeneracy, we find that, as energy is 
added to the system, p in equation (A9) remains essentially unchanged while T increases as c decreases. The input 
energy is mostly converted to an increase in temperature instead of doing the PdV work in a degenerate state. 
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Figure 10 depicts the coefficient of thermal expansion at constant pressure S 
(= —{d\np/dlnT)p), xt (= {din P/dlnT)p), and JS (related to the adiabatic index, see 
eq(Al8)) for a planet of 0.63Mj as a function of the co-moving mass coordinate m/Mp. The 
boundaries m/Mp = and m/Mp = 1 represent the location of the planet's center and its 
photosphere, respectively. The curves in Figure 10 are not smooth because the calculation involves 
several numerical derivatives. The magnitude of these thermal expansion coefficients decreases 
as the gas changes its phase from the ideal to the non-ideal regime, and approach to zero as a 
fluid increases its degeneracy, meaning that PdV work is of less importance in a more degenerate 
gas as indicated in equation (5). Figure 10 shows that 5 and XT monotonically increase with 
m/Mp, resulting from a transition from partial degeneracy in the inner region featuring the 
pressure-ionized hydrogen gas, to the non-ideal phase of dense molecular hydrogen in the middle 
range of m/Mp, and then to the ideal-gas regime in the very outer region (not able to be shown 
because of the small scale) where 5 = xt = 1- The magnitudes of 5 and xt are larger in the 
case of a relatively large size planet (with Rp = 3.22Rj) because the equation of state for its 
interior is less degenerate and more ideal than that for a smaller planet (with Rp = 1.66i?j). This 
correlation is consistent with the consequence that the temperature variation results from the shift 
of degeneracy. Since c drops as Rp increases, the term c^-^^ in equation (All) indicates that the 
exponent of Rp should become less than 10.32 as the degree of partial degeneracy is reduced. 



A. 3. Adiabatic interior composed of an ideal gas 

For planets with extremely large Rp{» 2Rj), the density near the center of the planet 

become sufficiently low that degeneracy is lifted and the equation of state is better approximated 
by that of an ideal gas for an n = 1.5 polytrope. In this limit, the radiative envelope is also 
relatively extensive. Integrating the radiative diffusion equation over the radiative envelope, we 
find the ratio of the temperature at the boundary Tb to the temperature at the photosphere Tph 
(Cox & Giuh 1968) to be 

n _[ {l + neff)Vph-l\ 



Tph \il +neff)V ad -'^ J 
where k oc p"^T~^, Ueff = (s -|- 3)/(m -|- 1), and the value of V at the photosphere 



(A12) 



V — ^ PphUph ^ _ ^ ^ ( 
'^^'^ - IGTracG T^, Mp ^^T^RlaT^^ 

In deriving the above equation, we have used Pph = 2GMp/3RpKph- Since C is usually much 
smaller than the stellar irradiation AirRpaTpi^ and therefore Vph << 1, the ratio T^/Tph is not a 
sensitive function of JC and Mp. Since Tph is determined by the stellar irradiation, neither it nor 
Tjj varies significantly with jC and Mp. 

There are some uncertainties for the values of m and s. If we evaluate k for the 
radiative/convective interface, we find from equation (A6) that m = and s = —4.08. The 
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temperature near the surface layer is 2000-3000 K such that the diatomic molecular hydrogen 
attains n = 5/2 and V^d = 2/7. With these parameters, we find from equation (A12), 
Tfy ~ 1.326Tp/i ~ 2150 K when Tp^ = 1620 K. Near the photosphere, the grains may also provide 
the dominant opacity source, in which case s = —1 but other parameters have the same values. In 
this case, we find from equation (A12) that l.91Tph ~ 3090 K which is too high for grain 

opacity to be relevant. The actual values of s and Tj, are probably between these extreme cases. 
Finally, with equations (A5) and (A6) we find for n = 5/2, an ideal-gas equation of state leads to 

C oc M^/^Rl/\ (A14) 



A. 4. Polytropic evolution and planet compressibility 

As we can see clearly from equations (A8), (All), and (A14), the values of 7 decrease from 
00 to 10.32 to 0.5 as the degeneracy is lifted for increasing values of Rp. The physical properties 
responsible for this change can be attributed to the "compressibility" of the planet. It is well 
known that the compressibility of a gas is linked to the effective exponent T^f f = dlnP/dln p 
which varies with different thermal conditions in a situation of thermal equilibrium. However, in 
the case of inflated giant planets, Fe// is always equal to the adiabatic value Fa^ = {din P/ din p)ad 
due to an efficient energy transport by convection. Consequently, Fe// varies with different 
equations of state which are characterized by the index of polytrope n in our case. We shall 
elucidate this point in this section. 

In general, equations (A5) and (A6) can be approximated with a simple relation 

Co,^. (A15) 

PbTb 

The factor Mp comes into the above expression because the larger Mp is, the larger gravity is 
in the thin radiative envelope, and therefore the stronger jC is due to the steeper temperature 
gradient in the radiative envelope. Unlike the relation with Mp, £. decreases as pb increases simply 
because the radiative energy flux decreases as the density and therefore the optical depth increase. 

The numerical solutions for the constant e models show that Ti, increases with Rp when the 
degree of partial degeneracy is relatively high. But Tb decreases with Rp when the degree of 
partial degeneracy is modest {Rp is larger than 2.3 Rj). This tendency suggests that Tb, a thermal 
quantity introduced by the equation of state in the framework of a polytropic analysis, primarily 

varies with c as shown in equation (A9). 

Figure 11 illustrates the mass density at the interface between the radiative envelope and 
the convective interior, logp^ as a function of the planet size log Rp. As the planet expands, 
log Pb decreases but its slope gets flattened. The similarity between Figure 11 and 1 suggests 
that the flattening of log pb is the major cause of the flattening of log £. as the planet is inflated. 
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Equation(A5) implies that 



l—n 




) 



( 



) 



n 



(A16) 



As demonstrated in Figure 4, the simulation shows that the fraction of the volume of the planet 
which deviates from the n = 1 polytrope increases with Rp. This dependence is roughly equivalent 
to the polytropic index being an increasing function of Rp. We fit the curves in the region of 
convection zones in Figure 4 for the core-less cases with a polytropic index and show the results in 
Table 1. The terms (7?.//x)" and in the above equation increase with n, but their influence is 
much weaker than the terms Rp~^ and Mp^~'^ because Rp and Mp are much larger. As a result, 
the simple polytropic approach roughly suggests that £ oc Rp^^"^ Mp^ ^ and hence the value of 7 
(as defined by eq. A8) decreases with Rp as — n decreases with Rp. Equation(A5) gives a more 
precise representation of 7 as a function of Rp-. 



The approximations made for the last expression in the above derivation are that we ignore the 
terms associated with the 2nd derivative with respect to log Rp, and that the terms c/log fe„/dlog Rp 
and dlogT^/dlogRp are less than unity. 

The magnitude of C increases with Mp for a given Rp as suggested by the factor Mp", which is 
consistent with the results shown in Figure 1. An increase of n with Rp in the polytropic analysis 
can be affirmed by the numerical results which show that the central density pc and the central 
pressure Pc do not scale as fast as Rp^ and Rp'^ respectively. But, pc oc Rp^'^ and Pc oc i?^^'^'^, 
i.e. the ratio of the central density to the average density pc/ < p > increases with Rp, meaning 
that n also increases with Rp. Hence, the fact that an inflating planet loses its partial degeneracy 
can be interpreted as an increase of n from 1 in the polytropic analysis. The overall effect on the 
intrinsic luminosity JC for its dependence on Rp is that 7 (= din jC / d\n Rp) decreases with Rp due 
to the flattening of the value of pbTj, as the planet inflates and loses its degeneracy. 

Another piece of information which suggests a correlation between electron degeneracy and 
7 is the comparison of the evolution of 7 between the core and core-less cases. As shown in §3, 
7 passes 5 at a smaller Rp for a less massive planet {Mp = 0.63, IMj) with a core than for that 
without a core. In the case of a less massive planet, a planet with a core requires more internal 
heating (see Fig. 1) than that with no core to be inflated to the same size, resulting in a lower 
degeneracy of the planet interior and therefore a faster decrease in 7 even though the core does 
not expand at all and the planet with a core is more gravitationally bound. On the other hand, 
7 decreases below 5 when Rp expands beyond 2.243i?j for a massive planet {Mp = 3Mj) in both 




-2 



dn 



(A17) 



d log Rp 
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core and core-less cases. The reason for this independence of the core structure is because a SMj 
planet without a core has a central density comparable to the 5 g/cm^ which we impose for the 
density of the core in our numerical prescription; i.e. the interior structure of a massive corc-lcss 
planet is comparable to that of the planet with a core, resulting in the similar C and D required 
for the core and core-less cases to be inflated to the same size (see Fig. 1). In summary, given a 
size of an inflated less massive giant planet and assuming that the core is not heated and does 
not radiate, the planet with a core is less degenerate than the one without a core, and hence the 
planet with a core has a smaller 7 and expands more for a given amount of entropy input (see 
equation (A21) and detailed explanations below). 

If we trace back through the above derivation, we note that the relation C oc R^~^Mp^ 
originates from equation (Al). The parameter K might be related to the total entropy content 
of a polytropic object. The energy equation can be expressed as follows (Kippenhahn &; Weigert 
1990): 

ds = CpJ (^d\nP- -^dlnpj , (A18) 

where J = — Va^ + 1/xt and xt = {dlnP/dlnT)p. Therefore the adiabatic exponent T and the 
polytropic index might be written in terms of S and xt- 

l + - = r = i-. (A19) 

n Jo 

For a completely degenerate gas, XT 0, 5 ^ 0, and 1/J5 can be written as the expression 
din P/ din p which equals 5/3 in the non-relativistic case. As the degree of degeneracy is reduced, 
both S and xt increase away from zero, and J, whose evolution is dominated by 1/xt, decreases 
from infinity. In the region of the middle and larger values of m/Mp where molecular hydrogen is 
abundant, both S and xt increase with Rp as the non-ideal effect is lifted (Saumon et al. 1995)^. 
This qualitative evolution of 6 and XT is in agreement with the plots obtained for a simulation 
with an inflated planet of 0.63 Mp (first two panels in Figure 10). The third panel in Figure 10 
shows that the evolution of the product J6 is dominated by 6 and therefore increases as the planet 
expands, leading to an increase in the polytropic index n. 

Motivated by the entropy equation (A18), which indicates that the specific entropy 
s oc ln{P/p^) = InK, we associate the issue of entropy to the polytropic approach with a constant 
n by re- writing equation (Al) as 

dlnRp = CMdlnK + CKdlnMp, (A20) 



^The effect that 5 increases as the molecular gas reduces its non-ideal properties can be seen from Fig 8 in Saumon 
et al. 1995 by noting that 6 = Xt/Xp: where Xp = {dlnP/dln p)t- The rise of xt above 1 during the inflation of 
a young giant planet is caused by the effects of molecular dissociation (see Fig 7 in Saumon et al. 1995 for more 
details). 
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where we have defined the compressibihties Cm at constant mass and Ck at constant entropy as 
fohows 



For a completely degenerate gas, K =constant (hence the result n/(3 — n) shown in eq. (A21) 
breaks down), Ck = — 1/3 (since n = 1.5), and the term associated with Cm vanishes in the 
non-relativistic case. In the cases of terrestrial planets and asteroids where the atomic/molecular 
interaction dominates and therefore the density is on the similar order regardless of mass, volume, 
and entropy, the interior structure can be roughly described by the n = polytrope. When n = 0, 
Ck = 1/3, Cm = and F = oo. These relations reflect a nearly constant density {Ck = 1/3) for 
the planetary interior as well as the difficulty of producing any non negligible density gradient with 
an entropy (Cm = 0) and gravity (F = oo) distribution. In the special case of partial degeneracy, 
Ck — and Cm increases as n rises from 1.^ The zero value of Ck indicates a maximum size of 
the planet at n = 1 solely in response to the mass change (Hubbard 1984). Phenomenologically 
this result arises from a transition from degeneracy {Mp oc R~^) to the non-ideal equation of state 
{Mp oc R^) due to atom-atom/molecule-molecule repulsion (Burrows et al. 2001; Shu 1992). In 
terms of the polytropic index, the state of partial degeneracy with n ~ 1 is just a transition from 
an n = 3/2 degenerate state to an n = constant-density state. The rise of Cm with n means 
that, for a given mass, the rate of Rp change in response to the change of total entropy gets more 
sensitive as n increases, leading to the radius-adiabat relation and therefore the radius-luminosity 
relation for constant heating per unit mass. 

It is physically straightforward to see why Cm increases with n and (therefore decreases 
with F). It is because the adiabatic exponent F = F^^ is just the bulk modulus for adiabatic 
expansion/compression. Equation (A21) can be re-arranged to have the following form: 



where we have defined the bulk modulus at constant planet's mass Fm^ = {d In Pc/d In pc)Mp 
which equals the constant 4/3. The adiabatic bulk modulus F^^ should decrease as the planet's 
elasticity Cm increases due to the reduction of electron degeneracy and non-ideal effects. Figure 12 

illustrates the isothermal curves for F^^ as a function of log p in the case of hydrogen (left panel) 
and helium (right panel). The number marked on each curve denotes the logarithmic value of 



^Equation(A21) gives the relation cICm /dlogRp = C\i{'i/n^)dn/d\ogRp > 0. A more precise 
expression derived from equation(Al) is given by the equation dCu /dlogRp = Cif [(6/n^)dn/dlog/Jp — 

{2/n^){dn/d\og Rpf \og{Mp / Rl.) + (l/n^){d^n/d(log Rpf)log{Mp / R%)+d'^ log kn/d{\og Rpf], which approximately 
equals C\i{Q/'n?)dn/d\ogRp when the other terms associated with d^ /dilogRp)^ and l/(log/?p)^ are neglected. 




(A21) 



(A22) 



^ = 3(F-FmJ, 



(A23) 
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temperature logr(K). The plots are drawn based on the tabulated equations of state (Saumon et 
al. 1995). The peaks of Tad (~ 2) around p ^ I g/cm^ at the temperature T ~ 10^-lO^K roughly 
correspond to the pressure-ionized regime of a Jupiter interior. This effect may be compared to 
the sudden rise of Tad for logT(K)= 3.22 at high densities as a result of non-ideal effects in the 
dense molecular (for hydrogen) or atomic (for helium) regime (also see Figs. 8 & 15 in Saumon et 
al. 1995). The variation of Xp = dlogP/dlog p\t and of Fad are quite similar; the rise of r^d in 
the pressure-ionized phase in the case of higher temperatures should also result from the non-ideal 
effects due to the interactions between densely-packed hydrogen and helium atoms, increasing the 
rigidity of the fluid. After all, equation (A20) describes the simple fact that the size of a planet 
is in general determined by its gravity (Mp), its entropy content (K), and the elastic properties 
of the planet responding to gravity (described by Ck) and to the entropy content (quantified by 
Cm)- 
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Table 1. Polytropic Fitting for An Inflated Giant Planet of IMj Without A Core 



Rp/Rj 


1.777 


1.973 


2.314 


2.711 


2.826 


n 


2.058 


2.141 


2.222 


2.358 


2.392 


(fen)" 


0.098 


0.102 


0.107 


0.118 


0.122 



Table 2. Parameters for various models: = 5 x 10~^Lq, Mp = IMj, and Tq = 1500 K. The 
last column entitled "figure" shows the figures, in terms of their labels, in which the model is 
plotted. For instance, Model 1 appears in Figures 5, 8, and 9. 
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log(Rp/Rj) log(Rp/Rj) 



Fig. 1. — Intrinsic luminosity as a function of planet's radius in logarithmic scales for three different 
masses. While the cases for planets without a core are displaced in the left panel, the cases for 
planets with a core are shown in the right panel. Discrete points arc the simulation data, which are 
connected by fitting curves associated with three different masses: 0.63Mj (solid line), IMj (dashed 
line), and 3Mj (dotted line). The number marked next to each discrete data point indicates the 
degeneracy D at the planet's center in the core-less cases (left panel) and the degeneracy D on the 
surface of the solid core in the core cases (right panel). All three curves get flattened as the planet's 
size Rp increases. 




P (Mbar) P (Mbar) 



Fig. 2. — The temperature profiles against the pressure inside a planet of IMj with (right panel) 
and without (left panel) a core. Different curves are plotted for different Rp: The interior structure 
evolves from the lowest curve to the upper-most one as the planet expands. 
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Fig. 3. — The mass density profiles as a function of the temperature for various planet's sizes in the 
case of having a core (right panel) and having no core (left panel) . The interior structure evolves 
from the upper-most curve to the lowest one as the planet expands. 




Fig. 4. — P/ (in c.g.s units) as a function of the co-moving radial coordinate m/Mp for various 
planet's sizes in the cases of a planet with a core (right panel) and without a core (left panel). 
m/Mp = and m/Mp = 1 represent the location of the planet's center and the bottom of its 

photosphere respectively. The interior structure evolves from the lowest curve to the upper-most 
curve as the planet expands to different sizes indicated by Rp/Rj shown at the upper-left corner 
of each panel. 
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Fig. 5. — The density profiles of a giant planet in the cases of Model 1 &: 6 (dotted curve), Model 2 

(solid curve), and Model 5 (dashed curve). Model 1 is a case of uniform heating, while model 6 is a 
case with the same integrated heating rate but with the energy deposition concentrated in a shell 
at mass fraction mo/Mp = 0.05. Models 2 and 5 also have the same heating rate but the energy is 
deposited in shells of mass fraction mo/Mp = 0.9999 and 0.70, respectively. 
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Fig. 6. — The evolution of Rp, the temperature profile, and the density profile for Models 11 (3 
upper panels) and 12 (3 lower panels) in Table 2. While we only show the data after the planet 
has reached a quasi-thermal equilibrium in Model 11, in Model 12 the planet is not in the quasi- 
thermal equilibrium until the data point 4 (see the lower- left panel). The three solid curves from 
right to left in each of the temperature and density plots for Model 11 (Model 12) correspond to 
the three different stages marked by 1, 6, and 11 (1, 4, and 10) respectively in the Rp vs t plot. 
The vertical dashed lines in the temperature and density plots mark the location of the maximal 
heating of the Gaussian heating profile for Model 11 {mo/Mp = 0.9) and Model 12 {mo/Mp = 0.7). 
The temperature and density profiles for the uniform heating per unit mass represented by Model 
9 (dash-dotted line) and for the surface heating denoted by Model 10 (dotted line) are also plotted 
for comparison. 




Fig. 7. — The radial profiles of J 5 (= l/F) for the data point 1 (solid curve) and the data point 
4 (dotted curve) in Model 12. The data points 1 and 4 are indicated on the lower-left panel in 
Figure 6. In Model 12, the energy is deposited in the shell of mass fraction mo/Mp = 0.7. 
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Fig. 8. — The comparison of the density profiles for the regular grain opacity na = kq (the solid 

curve for Model 1 and the dashed curve for Model 4) and for the reduced grain opacity = 10~^ko 
(the solid curve for Model 13 and the dashed curve for Model 14). While the heating is deposited 
uniformly in mass in Models 1 and 13, Model 4 and 14 have the same integrated heating rate but 
with the energy deposition concentrated in a shell at mass fraction nio/Mp = 0.9. 
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Fig. 9. — Temperature and opacity profiles in the radiative envelope for Model 1 (solid curve), 
Model 16 (dashed curve), and Model 17 (dotted curve). Different values of irradiation temperature 
are imposed here: Tq = 1500 K for Model 1, O.STq for Model 16, and O.bTo for Model 17. 
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Fig. 10. — The coefficient of tliermal expansion 5, xt, and JS as a function of the co-moving mass 
coordinate m/Mp, plotted here for an inflated planet of 0.63Mj with two different radii: 1.66-Rj 
(solid line) and 3.22Rj (dotted line). 
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Fig. 11. — The density at the interface between the radiative envelope and the convective interior 
Ph as a function of the planet's radius Rp. While discrete points are the simulation data, the fitting 
curves are illustrated by a solid line (0.63Mj), a dashed line (IMj), and a dotted line (3Mj). 




Fig. 12. — Tad as a function of log/?(g/cm^) and logT(K) for hydrogen (left panel) and for helium 
(right panel). The number marked on each curve denotes logT(K). 



